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Abstract. Tribological aspects must be taken into account for a sustainable design 
of new components and materials developed to obtain weight reduction and greater 
efficiency. Reducing friction and wear produces energy and material savings, both 
connected with several Sustainable Development Goals. To limit the time consum- 
ing expensive experimental tests on new materials and components, simulations 
can be performed for which reliable values of the friction coefficient are necessary. 
In this work, some basic aspects of the lubrication regimes are firstly reviewed with 
the related friction coefficient trends represented with the Stribeck and Lambda 
curves, also evidencing the reasons of the similarity between the two curves. For- 
mulas and diagrams are then reported for the friction coefficient of full lubricated 
conformal pairs. For thrust bearings the friction coefficient f can be expressed 
as a function of the parameter m and is related to the Kingsbury number K. For 
tilting pads f is proportional to K the power of 0.5. For plain journal bearings f is 
a function of the dimensionless eccentricity £ and is related to the the Sommer- 
feld number S to powers ranging roughly from 0.5 to 0.8 depending on the ratio 
between the axial length and the diameter for § smaller than 0.1, and tending to 1 
for higher values of S. The reported formulas and diagrams can be used for design 
purpose. 
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1 Introduction 


Tribological aspects must be taken into account for a correct design of new components 
and materials. Reductions of friction and wear are connected to greater efficiency and 
material saving. This means less pollution and therefore a more sustainable life, strictly 
connected with the Sustainable Development Goals, where efficiency is mentioned sev- 
eral times [1]. Particularly Goal 12 (Responsible Consumption and Production) mentions 
sustainable consumption and production, efficient use of natural resources, environmen- 
tally sound management of chemicals and all wastes, and the significant reduction of 
release to air, water and soil. There is a trend today to weight reduction and more 
power density of the mechanical components. New materials, surface treatments, coat- 
ings and textures can be used and their tribological behavior under lubricated conditions 
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must be investigated. Experimental tests are necessary due to the complexity of the 
tribological phenomena, but they can be extremely energy- and time-consuming. Theo- 
retical/numerical simulations of the mechanical components are very useful for reducing 
the number of the expensive experimental tests. A reliable value of the friction coefficient 
is necessary for an appropriate simulation of the behavior of the pairs. Unfortunately 
there are a lot of factors that influence friction, as the lubrication regime, the behavior 
of the lubricant and the geometry of the contact (e.g. conformal or non-conformal). Due 
to the vastness of the problem, this work is necessarily focused only on some aspects, 
in particular on friction in lubricated conformal contacts. After a recall to some basic 
aspects, friction coefficient formulas and diagrams under full lubricated conditions for 
conformal pairs with Newtonian lubricants are reported that can be used for a sustainable 
tribological design. 


2 Friction Coefficient Trends: Stribeck and Lambda Curves 


The values of the friction coefficient in lubricated contacts can vary a lot depending on 
the level of occurrence of contacts between the roughness asperities of the bodies’ sur- 
faces. Three different regimes are commonly identified: boundary lubrication, when the 
lubricant film is reduced to molecules absorbed on the solids’ surfaces, mixed lubrica- 
tion, when the load is shared between the locally contacting surfaces and the pressurized 
lubricant film, and full fluid lubrication when the lubricant film completely separates 
the solids. The transition among the different regimes is well depicted by two similar 
curves, the Stribeck and Lambda curves, Fig. 1. 


Fig. 1. Qualitative trend of the friction coefficient in Stribeck and Lambda curves. 


As reported in [2], in the Stribeck curve the friction coefficient is usually plotted as a 
function of a service dimensionless parameter Sp proportional to the lubricant viscosity 
u and to the body’s surface velocity u, and inversely proportional to the load per unit 
length W/L, Sp x (uu)/(W/L). This kind of friction trend was firstly experimentally 
observed in journal bearings by Stribeck in 1902, as reported in [3]. Different expressions 
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and names can be found in literature for Sp, mainly depending on the kind of lubricated 
contact [2]. Different curves can be obtained by varying the surface roughness. The 
roughness effect is instead included in the Lambda (A) curve, where the friction coef- 
ficient is plotted as a function of the ratio of a representative value of lubricant film 
thickness and the equivalent root mean square roughness R, of the surfaces of the two 


bodies in contact, A = h/ R? Lt Rio: There is a connection between Sp and A and 


this explains the similarity between the two curves as reported below. 

A comprehensive expression of the friction coefficient f can be found by using 
the ratio between the tangential (friction) and the normal force acting between two 
contacting bodies including both the effect of the solids’ contacts and of the viscous 
action of the lubricant. By indicating with fp the friction coefficient under boundary 
lubrication conditions and with f, the friction coefficient for full lubrication conditions, 
the total friction force can be expressed as T = fyWp + fn Wh, being W = Wp + Wp the 
total load, which leads to the total friction coefficient formula [3]: 

T Wp Wh 
f= z h thy (1) 
with W/W and W;/W ranging from 0 to 1. When the ratio W/W = 1 (Wa/W = 0) 
boundary lubrication conditions occur and when W;/W = 0 (W/W = 1) full lubrication 
conditions occur. Different expressions can be found in literature for the ratio W;,/W as a 
function of the dimensionless film thickness A. By indicating this load sharing function 
as g(A), Eq. (1) can be put in forms such as 


f =f[1 —e(A)] +h g(a) (2) 


The values of g(A) range from 0 for A = O up to 1 when A reaches the value 
corresponding to the boundary between mixed and full lubrication conditions; g (A) 
maintains the value of 1 for A > 1. Expressions for g(A) can be found for instance in 
[4]. 

Regarding with the values of the friction coefficients, fp is usually considered con- 
stant (typical values range from 0.08 to 0.1), while f} has to be evaluated from the friction 
force generated by the fluid action on the body’s surface. Actually things are more com- 
plicated: in boundary and mixed lubrication regimes wear can occur that affects fp; the 
effect of roughness on the film thickness under mixed lubrication conditions depends 
on additional aspects rather than R4 only, as the orientation of the machining marks; 
the calculated lubricated friction depends on the lubricant’s constitutive behavior that 
for non-conformal contacts is often non Newtonian, depending on the contact geometry 
and on the working conditions. Every aspect needs specific studies. The results reported 
in the following are specifically addressed to friction for conformal contacts under full 
fluid lubrication conditions with Newtonian fluids. In this case the friction force is the 
integral of the tangential stress Tt: 


where u is the dynamic viscosity, v the velocity of the fluid and y the shear strain rate. 
As it can be also noted from Fig. 2, dv/dy is variable with x (in addition than along y). 


202 E. Ciulli 


To estimate the action on a body, t must be evaluated at the surface level (e.g. y = 0 for 
the moving body in the figure). 


Fig. 2. Generic lubricated contact. 


A rough estimation of the friction coefficient f p can be obtained as the ratio between 
approximated mean values of shear stress and pressure p: 


fa = — (4) 


where t* is the value of the tangential stress in the position x* where the pressure is 
maximum and the film thickness is h* (Fig. 2): 


u 


T* = M (5) 


and pm is the mean contact pressure. For the full lubricated contacts, the film thick- 
ness is always a function of the viscosity, the body’s velocity, a significant geometrical 
dimension L and the load: 
peo 1 uo2 LS 
x ————————s 
wes 


h* (6) 

The values of the exponents c1, c2, c3 and c4 depend on several factors as the kind 
of lubrication and the shape of the contacting bodies. The mean pressure is related to 
the normal force W and to the pair’s dimensions (let’s indicate them with a generic R): 


ws 
Pm X ee (7) 
By combining Eqs. (4), (5), (6) and (7), we obtain: 
t* puw’s R“ ulule? Ro 
= oe we =a le (8) 
Pm uetus? L3 Wes (W /L)°S~ © | Le3+es—c4 


The necessary constant that must appear in an equality is generally dimensional. 
The expression in the square brackets contains the same quantities present in Sp and in 
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the film thickness formulas, but with different exponents. This explains the similarity 
between Stribeck and Lambda curves. In fact, under full lubrication conditions we have: 


Ci yë? 1 
os ae / -a 

2 2 

Ra + Ro 
If for instance c1, c2, c3 and c4 are equal to 0.5, cs = 1 and cg = 2, as for the 


thrust bearings described in the next section, we obtain (the necessary constants for the 
equalities have the dimensions of an inverse of a length and of a length respectively): 


0.5,,0.5 R2 R2 0.5,,0.5 1 K 
maa E a vo Aa IS —— « ` È = (10) 
(W/L) WD [R + Re, JR +R, 


3 Thrust Bearings 


Expressions for W and T acting on a single rectangular pad (a parallel-surface slider 
bearing, Fig. 3) as a function of the significant quantities of the lubricated contact can 
be found in several papers and books. The friction coefficient f p can then be evaluated 
by dividing T and W. Analytical expressions can be found for the different functions 
when b > a, while they can be numerically evaluated in the other cases. For the sake of 
brevity, only the relevant formulas for the evaluation of f p are reported in the following, 
while for instance [5] and [6] can be seen for the additional expressions. 
The Kingsbury number K is commonly used as Sp for the pads of thrust bearings: 


b 
SO alsa laces (11) 
Pma W 
with pm = W/ab the mean pressure acting on the pad. 
By using the dimensionless quantities (refer to Fig. 3 for the symbols used): 
hi — ho 
m= , A= b/a (12) 
No 
the following expression can be found for the friction coefficient: 
T 
fr = = VKo(m, A) (13) 


For a given tilting pad, à can be calculated from the dimensions a and b, while m is 
fixed. Therefore f p depends on VK. Expressions for the function ¢ are reported below. 


3.1 Infinitely Wide Slider Bearing 


In the case of b > a the side-leakage along the z direction can be neglected and the 
flux is one-dimensional. Analytical expressions can be found for the main quantities as 
a function of the parameter m only. Rearranging the formulas reported for instance in 
[5], the friction coefficient function can be expressed as: 


2 (4+ 2m)ln(1 + m) — 3m 

ọm) = y (14) 
3 /(2+m)[2 + m)In( +m) — 2m] 

The trend of this function is shown with a dashed line in Fig. 4. 
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Fig. 3. Parallel-surface slider bearing. 


3.2 Finite Width Slider Bearing 


Curve-fitted functions based on numerically calculated results are available in [7] for the 
quantities of interest for the design of hydrodynamic bearings. The Reynolds equations 
is numerically solved for certain values of the input parameters, so that the expressions 
found for the rectangular tilting pad thrust bearing are valid in the domains 0.8 < m <5 
and 0.5 < à < 6. By rearranging the formulas of T and W the following expression can 
be obtained for the friction coefficient function: 


__ fom% > 


15 
1+2m (a) 


m, A)={1 
$m, A) ( a 
where fọ = 0.7095 + 0.2873e7^, gg = 1.1964 + 0.0715e7^, fy = 1.621 — 
0.5e719964, gy, = 0.4607 — z7% and ky = 0.6709 + 20105 , 

Some trends of ¢ (m, à) are shown in Fig. 4. It’s worth noting that values for b/a = 
4 are very close to those obtained with Eq. (14) for the infinitely wide slider bearing. 


Fig. 4. Trends of the friction function ¢ for slider bearings. 
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The Sommerfeld number S is normally used as Sp for plain journal bearings, Fig. 5: 


ç- Q (DY) _ w2LD/ DY" it 
~ 20pm (2) ~ 2nW (2) oe) 
where {2 is the rotational speed of the shaft, D and L the bearing’s diameter and axial 
length respectively, pm the conventional mean pressure W/(DL), C the radial clearance. 
It is worth noting that, as reported in [8], the original dimensionless Sommerfeld number 
was in reality the reciprocal of the one commonly used shown in Eq. (16). 


Results are usually presented as a function of the dimensionless eccentricity (being 
the eccentricity e the distance between the centers of shaft and housing): 


e=e/C (17) 


The friction coefficient is evaluated as the ratio of the friction torque My and the 
product between the load and the bearing’s radius: 


Mr 


TIP (18) 

Both S and fp can be evaluated as a function of the dimensionless eccentricity 
£. Expressions are available for long (with different boundary conditions) and short 
bearings, while tabulated numerical results are available for intermediate values of the 
ratio L/D. In fact, when the dimension (in this case the axial length L) along the direction 
perpendicular to the one of the motion is much bigger or much smaller than the other 


dimension, an analytical solution can be found. 


Fig. 5. Plain journal bearing. 


4.1 Long Journal Bearing 


The solution for long bearings, usually valid for L/D > 4 [6], is found when the side leak- 
age along the bearing’s axis direction can be neglected. Different boundary conditions 
can be used that produce different expressions for the main quantities. The expressions 
reported below were obtained by rearranging the ones reported in [6]. 
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Sommerfeld’s Boundary Conditions. The so called Sommerfeld’s boundary condi- 
tions only consider the periodicity conditions that the pressure p is the same for @ 
= 0 and g = 27, where the film thickness is maximum (see Fig. 5). This produce an 
anti-symmetrical pressure distribution respect to the position with g = x. If p(0) = p(27r) 
= 0, p < O for 2 < < 2 x. The following expressions are found for both S and f p: 


(2+6*)V1—e? | 201428? 
ne "DOB 


= (19) 
Giimbel’s or Half Sommerfeld’s Boundary Conditions. By using the Sommerfeld’s 
boundary conditions and simply putting to zero the pressure for 7 < p < 27, we obtain: 


(2 + 7)(1 —«7) 2C (4+ 2e + Se? + 463) V1 — e? 
S= Sh = E (20) 
one 4e? + n?(1 -—£?) 6e(1 + e) 4e? + 22(1 — £?) 


4.2 Short Journal Bearing 


Analytical expressions can be also found when the axial length of the bearing is much 
smaller compared to its diameter (for L/D < 1/8 according to [6]). This is the so called 
Ocvirk and Dubois solution. The half Sommerfeld’s boundary conditions are also 
applied. By rearranging the formulas reported in [6], it is found: 


(PJ po (oje 
= s Jh = 
L ne fn? - e? (n? — 16) DAL eja? — e? (n? — 16) 


(21) 


It’s worth noting that in this case f , is not only related to 2C/D as for the long 
bearing but also to the ratio (DIL)*. 


4.3 Finite Length Journal Bearing 


For 0.25 < L/D < 4 the numerical solution of the complete Reynold’s equation is 
necessary. Tabulated numerical results for the main design quantities can be found in 
[6] for L/D = 0.25, 0.5, 1 and 2 for some values of £ between 0.1 and 0.95. Values are 
reported for the friction coefficient function: 


‘= p 22 
f = faz (22) 


Values for L/D < 1/8 and L/D > 4 are also reported in the table, corresponding to 
the short and long bearing solutions respectively. In the first case S and f’ are reported 
multiplied by (L/D) in order to obtain a function of ¢ only. In the second case the same 
results can be obtained with the solution for long bearings by imposing the so called 
Reynold’s boundary conditions p(0) = p(~o) = 0 where the angle go can be evaluated 
with the additional condition dp/d¢|yg—g) = 0. 
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Fig. 6. Trends of the friction function f’ for plain journal bearings. 


In Fig. 6 f” is plotted a function of S for the different L/D ratios together with the 
results for long and short bearings. Values of f’ increase by decreasing L/D with an 
increasing value of the exponent of S. For a value of S lower than 0.1 f’ « S”, with n 
depending on L/D and roughly ranging from 0.5 for the long bearing solution up to 0.8 
for the short one. For greater values of S the trend of f’ tends to become linear. It can 
also be noted from its irregular trend for S < 0.03 that the full Sommerfeld solution can 
be used only for very light loads (high values of S). 


5 Conclusions 


The friction coefficient trend is well represented by the Stribeck and Lambda curves. In 
both diagrams the friction coefficient is plotted as a function of dimensionless quantities 
containing viscosity, velocity and load, and this explain the similarity between the two 
curves. Formulas and diagrams are reported for the friction coefficient of full lubricated 
conformal pairs with Newtonian fluids. Based on analytical and numerical solutions, 
trends of the friction coefficient for thrust and plain journal bearings are reported. For 
tilting pad thrust bearings the friction coefficient is a function of the parameter m and 
is related to the Kingsbury number to the power of 0.5. For plain journal bearings the 
friction coefficient is a function of the dimensionless eccentricity € and is related to the 
the Sommerfeld number to powers ranging roughly from 0.5 to 0.8 depending on the 
ratio between the axial length and the diameter, and tending to | for high values of S. 
The formulas and diagrams can be used for a quick estimation of the friction coefficient 
for design purposes. 
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